Level statistics of quantum dots coupled to reservoirs 
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We study the effect of electron tunneling on the level statis- 
tics of quantum dots. While the coupling between individual 
levels and the electron reservoir leads predominantly to the 
expected level broadening, the indirect coupling of adjacent 
levels via the reservoir results in a new asymptotic level statis- 
tics and broadening. These features, which are attributed to 
renormalized spectral weights rather than renormalized eigen- 
values of the Hamiltonian, are observable in the asymptotic 
frequency dependence in an absorption experiment. 

Introduction. The statistics of the excitation spectra 
determines many physical properties of quantum dots. 
These include thermodynamic properties, transport co- 
efficients, the response to external a.c. fields, and statis- 
tical fluctuations For a model of non- interacting 
electrons confined within a moderately disordered quan- 
tum dot (diffusive disorder), the description of the single- 
particle levels for energies below the Thouless energy con- 
forms to the Wigner-Dyson theory of random matrices 
(RMT) §|-§. As a result of the level repulsion, on en- 
ergy scales smaller than the average single-particle spac- 
ing A, the statistics of the inter level spacings for s <C 1 
is given by P {s) = tts/2, if the random system pos- 
sesses time reversal symmetry, (orthogonal case) and by 
P u (s) — 32s 2 /tt 2 otherwise (unitary case) Here s is 
measured in units of A. This level statistics manifests it- 
self in the response to low-frequency fields. For instance, 
the a.c. dissipative conductance scales as the applied fre- 
quency wo, or as w 2 , for the two symmetry classes, respec- 
tively. This result, valid for the grandcanonical ensemble, 
is modified for the canonical ensemble Jl|,p|,p| JTc| , pd| ] . 

If the quantum dot with discrete bare levels is cou- 
pled to external degrees of freedom, these levels acquire 
a finite width T. While the broadening is known to de- 
pend on the coupling and on the energy measured from 
6f (see e.g. |l^]) it is commonly assumed that the lev- 
els statistics is unchanged. This means, the standard 
Kubo-Greenwood formalism can be employed, provided 
that the (5-function describing the energy conservation is 
replaced by a Lorentzian with finite width. 

In this Letter we study the influence of a dot-reservoir 
coupling on the effective level statistics and the asymp- 
totic wo-dependence of the absorbed power. A simple, 
though generic, model, where the quantum dot is cou- 
pled to an ideal reservoir of non-interacting electrons is 
used. In most of our analysis we ignore interactions in 
the dot. Breaking the RMT invariance by the choice of a 
preferred basis may modify the level statistics as it was 
demonstrated for disorder in metals in Rcf. Ol. In our 



case, the system is not invariant under rotation between 
the dot and the lead states. We show that the coupling to 
the lead gives rise to a renormalization of both the low- 
frequency spectrum and the bare level broadening T. We 
attribute the origin of the modified statistics and broad- 
ening (Eqs. to renormalized spectral weights 
rather than renormalized eigenvalues of the Hamiltonian. 
Subsequently the frequency- and temperature-dependent 
absorption is calculated, (Eqs. @ with @-Q). Fi- 
nally we comment briefly on the effect of Coulomb inter- 
actions. 

Model. The model Hamiltonian 
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describes a quantum dot with two levels a — 1, 2 coupled 
to the same electron reservoir fl4j| . The tunnel matrix 
elements are assumed to be independent of the reservoir 
state, Tfccr = T a . The Coulomb interaction will be dis- 
cussed at the end this Letter. 

Due to the coupling to the reservoirs, if both dot levels 
are well separated, each level acquires a finite width re- 
lated to r (T (w) = 2nJ2k \Tka\ 2 5{to — e k ). However, if the 
levels are degenerate, e\ = e 2 = e, a rotation within the 
dot subspace yields 



h' = E <*4**+eE c '° c ° + E [ T ^A 



(2) 



where c[ = {T lCl +T 2 c 2 )/ ^{T^ 2 + |T 2 | 2 and ^ = {T 2 c x - 
T1C2) I \J '|21| 2 + I X2 1 2 . In this form one recognizes that 
only the new level 1 is coupled to a reservoir of electrons 
(leading to a Lorentzian form of the spectral density) 
while the new level 2 is isolated (showing therefore a 5- 
peaked spectral density). Therefore, in the unrotated 
basis c\,c 2 , the spectral density of each level consists of 
a Lorentzian part but also contains a (5-function. The 
sharp peak persists if the degeneracy is slightly lifted. 

Spectral density. In order to catch the evolution of 
the sharp peak in the spectral density quantitatively, we 
employ the equations of motion for the Green's functions 
G a a(t) ■— ~iQ(t)({c a (t),c a (0)}). Since the Hamiltonian 
Eq. (nl) is quadratic, we can find the exact solution 



Efjfj En-r 



(3) 



with S^/ = J2 k T ka T^, /(u> - e k ), i.e. Y, aa (u + i0+) = 
-iT a /2 and E CTff {to + iO+)Y. Sa {uj + i0+) = -r CT r ff /4 (for 
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simplicity we put U = 1). The spectral function follows 
from A a (u) = — lmG aa (u) + iO + )/ir. For a ^ e 2 it is 



T a (u - e ff ) 2 /2tt 
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e ff ) 2 + [(r 1 + r 2 v/4- e f 



(4) 



with e = (rie 2 + r 2 ei)/(ri + T 2 ). In Fig. § we plot the 
spectral density of level 2 for several ratios T/Ae. For 
the following we define Ae = e 2 — t\. If the bare levels 
are well separated, i.e., Ae 3> T, the spectral density is 
approximately a Lorentzian with width r CT , 



A a {uj) = 



1 



r CT /2 



(5) 



and the maximum is at the bare value e a , i.e., the level 
splitting remains unrenormalized Ae = Ae. We see, how- 
ever, from Fig. [j] that things are different in the opposite 
limit Ae <C T. Then, the level position is renormalizcd 
which leads to a modified level statistics, and the width 
is much smaller than the bare broadening. 
For the degenerate case ei — e 2 we find 



r ff /27r 



(^-e CT r + (r 1 + r 2 )2/4 



and recover the expected S peak. 

Modified level statistics. For Ae <C T, Eq. 
tively reduces to the sum of two Lorentzians, 



(6) 
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It shows a broad peak with width Ti + T 2 , and a sharper 
one with width f (see inset of Fig. [I]) . For the latter peak 
the level splitting Ae and the width T are given by 



Ae 
Ae 
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(8) 



This result implies a new statistics for the renormalized 
interlevel spacing s — Ae/A in the regime s <C T/A <C 1, 
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(9) 
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and a renormalized broadening T — 4Fir 2 (Ae) 2 /(ri + 
r 2 ) 3 . In the orthogonal case the distribution even di- 
verges for s — > with integrable divergence. 

Results for a wider range of energies are displayed in 
Figs. ^| and ^ for the renormalized level spacing as a func- 
tion of the bare one and for the new level statistics, re- 
spectively, and compared to the low-energy asymptotic 
results (|), (|) and ©. 



What is the origin of the modified -P(s)? If we omit the 
tunneling part of the Hamiltonian Eq. ([j]) the the eigen- 
values, efe and Co-, have spectral weights totally within the 
corresponding states k and a. It is tempting to assume 
that the coupling of the level a to the reservoir states k 
leads to a renormalization of the eigenvalues e CT , thus ex- 
plaining the modified level statistics. But this is not the 
case. The crucial point is rather the renormalization of 
the eigenvectors. Due to the tunneling each eigenvector 
acquires a finite overlap with the bare dot states. What 
we see in Fig. |l| is the envelope of the square of this over- 
lap. The relevant physics is, therefore, a spectral weight 
renormalization effect rather than energy renormaliza- 
tion. 

Signatures in an absorption experiment. How can the 
evolution of the sharp peak be probed? In the following 
we discuss consequences for an absorption experiment. 
For a single quantum dot the absorption power of a pho- 
ton with energy u>q, accompanied with a transition from 
level a to level a', is given by 



da; A a (us - y) A <y' ( w + y) x 

-M)] ■ (ID 



Since in realistic experiments an ensemble of dots is 
probed, the total signal I aa i is determined by the average 
over all possible configurations, i.e., we have to average 
over the chemical potential f-^&n % ■ ■ ■ ~ ^ ■ ■ • 

as well as over bare level separation ds P(s) .... After 
performing the integral over d\i we get 



lav' 
D aa > = 



CJ 



A e^<> - 1 
1 

ds P{s)B, 



dioA a (oj - y) A <y' ( w + 



(12) 

(13) 
(14) 



The quantity B aa ' probes the properties of the transition 
for a given level spacing. It accounts for the renormal- 
ization of the level spacing and width. The total signal, 
however, is determined by the average over all (bare) level 
spacings, D aa i . For a = a' the absorption is accompa- 
nied by a transition within one level and for a =/= a' a 
transition between the two levels. In the following we 
are interested in the asymptotic behavior for oj — > 0, 
i.e., ujq < T < A. 

In order to establish some reference frame we first ap- 
proximate A a (oj) by (i) delta functions (which neglects 
the coupling of each level to the reservoir) and by (ii) 
Lorentzians with width r CT (which neglects the indirect 
coupling of the levels) before (iii) we use the exact form 
Eq. @. 
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(i) In the first case we find D\2 — P(ujq), i.e., I? 12 
7T 2 w /2 and Df 2 = 32^/vr 2 , while D xl = D 22 = D 2 i 
0. I.e. the absorption power vanishes for luq —* 0. 

(ii) In the second case 



B\2 — 



B 2 i = 



1 



2tt (cj - Ae) 2 + (Ti + r 2 ) 2 /4 

J_ r! + r 2 

2tt (w + Ae) 2 + (ri + T 2 ) 2 /4 

1 r„ 



2tt 0Jq 



r 2 ' 



which yields for luq = 



D o - D o _i r 1 + r 2 2A 
^i2-^2i- 4 A2 ln ri + r2 

^12-^21 - -3 A2 



(15) 
(16) 
(17) 

(18) 
(19) 



and D° a = 1/4IV, = 32/37r 3 r (T . Two points are 

remarkable. First, due to the broadening there is always 
an overlap of the shifted spectral functions, i.e. even for 
luq = 0, D±2 remains finite. Second, for systems with time 
reversal symmetry the (Ae) -2 behavior of Eq. ( [l5"| ) leads 
to logarithmic behavior, which is cut off at low energies 
by the average level width (Ti + T 2 )/2. 

(iii) The renormalization of the level splitting and 
width due to indirect coupling of the levels to each other 
is fully included in the exact spectral function Eq. (Ji|). 
The leading terms are 



(ri+r 2 

41^2 



u \2 — u 2\ 



8 ri + i 2 
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and 



4r a 8 a 2 r„ 



(ri+r 2 ) 3 
4Tir 2 



(20) 
(21) 

(22) 
(23) 

We note that the low-energy cutoff \J uj a is fre- 

quency dependent and determines the asymptotic ab- 
sorption power as can be seen in Fig. ^. In comparison 
to Eqs. (18) and ( |Tg| ) a factor 1/2 arises due to details in 
the spectral density. 

We remark here that for a symplectic ensemble we get 
similar result as for the unitary case Eq. ( pl| ) and (23) 



32 i6 r 1 + r 2 r g 
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21 -(2 16 /3V)(r 1+ r 2 )/A 2 a nd^ 

2 17 /3 7 7r 4 )(r ff /r CT )(r 1 + r 2 )/A 2 . 

Interaction. To get a qualitative understanding for the 
effect on interactions we include a finite charging energy 
in the Hamiltonian H — > H + XJn\n2- As an example 
we consider the degenerate case Ae = 0. A rotation 



within the dot subspace, as discussed above, can still 
decouple one new level from the reservoir. But the lev- 
els still influence each other due to the charging energy 
term U n^n^ and the conclusion that there will be a delta- 
function peak in the spectral density is no longer valid. 
If we neglect all terms in the equations of motion which 
would correspond to correlations in the reservoirs we get 
a closed set of equations for the Green's functions. The 
resulting expression for G aa (u> + i0 + ) for arbitrary U is 
too lengthy to be presented here. However, for U — > 00 
it simplifies to 
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(n) 



1 



lu — e — 2Ti(u>) 



(24) 



with ImE(w + i0+) = -T/2f(u) and ReE(w + i0+) = 
r/27r[ln(/3C//27r) - Re*(l/2 + i/L/27r)]. Here, $(z) is 
the digamma function and f(co) is the Fermi function. 
If the levels are far (in units of r^) from the Fermi 
level they only contribute to the absorption power in 
the case \u> — e\ < fcsT since otherwise the levels are 
totally empty or filled. Then, the second term acquires 
a width of the order of T and does not show a sharp 
peak. If the levels are close to the Fermi level also the 
low-temperature regime is important. But then the spec- 
tral density shows a richer structure related to collec- 
tive many-particle states. Similar effects lead to Kondo 
physics in a single-level spin-degenerate quantum dots 
with spin conservation (see e.g. |jl],[l(|). The analysis of 
the level's properties and their effect on the absorption 
power becomes more complicated and is out of the aim 
of this Letter. 

In summary, our model analysis produced a modified 
statistics of the broadened levels due to a renormaliza- 
tion of the spectral weights. As a consequence, in the or- 
thogonal case the asymptotic absorption power depends 
logarithmically on the external frequency. 
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FIG. 1. Spectral density for level 2 with Ti = T 2 = T and 
ei = —62 = 1. For Kea peak with width r < T evolves 
near zero. The spectral density of level 1 is related to that 
of level 2 by Ai(lj) = Az(— u). Inset: Spectral density for 
T/e2 = 25 (solid line) showing a sharp peak sitting atop a 
Lorentzian with width 2P (dashed line). 
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FIG. 2. Renormalized level spacing as a function of the 
bare one in double logarithmic plot for Ti = T2 = T. The 
dotted and the dashed curve are the high- and low-energy 
(see Eq. (^)) asymptotes, respectively. 




FIG. 3. Solid lines: new level statistics P(s) in the or- 
thogonal (upper curve) and unitary ensemble (lower curve) 
with ri = T2 = r = A/20. Dashed lines: low-energy asymp- 
totes Eqs. (^|) and ([lo|). Dotted lines: high-energy asymptotes 
(Breit-Wigner distribution) . 
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FIG. 4. For the orthogonal ensemble the asymptotic ab- 
sorption power is logarithmically divergent at small frequen- 
cies (solid line). In the simplified model where the indirect 
coupling between the levels is neglected, the low-energy be- 
havior shows a saturation at F (dashed line). The dotted 
curve is the asymptote Eq. (po|). 



